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Abstract 

We develop a method for the calculation of vacuum expectation values of local 
operators in the Lieb-Liniger model. This method is based on a set of new identities 
obtained using integrability and effective theory ("bosonization") description. We use 
this method to get an explicit expression for the three-body local correlation function, 
measured in a recent experiment 

1 Introduction 

In this paper we present an approach to the calculation of the vacuum expectation 
values of local operators in the Lieb-Liniger model. Although the technique is general 
we shall concentrate on an particular operator, {x)^ip{x)^ , whose vacuum expectation 
value was measured in a recent experiment H]. The physical motivation for the study 
of this problem and the analysis of the resulting expression has been given in Ref. Q. 
Here we focus on the mathematical aspects of the problem and present the details of 
the derivation. In some integrable models used in condensed matter physics the exact 
expressions for expectation values of local operators are known, Refs. ^. However, 
it is not known how to use these results for the calculations in the Lieb-Liniger model. 

The Lieb-Liniger model describes a one-dimensional gas of bosons interacting via 
a (5-potential ||5|. Its Hamiltonian is 



L 

H = I dx 







-^\x)d'^'4){x) + cil)"^ {x)iIj\x)'4}{x)'4}{x) , (1-1) 



where the units have been chosen such that h = 2m = 1. The boson fields ip and V'^ 
satisfy canonical equal-time commutation relations: 

[il^{x),tlj\y)]=5{x-y), [,p\x),tl;\y)\ = [^{x),i^{y)] = (L2) 



*Department of Physics, Lancaster University, Lancaster, LAl 4YB, United Kingdom 
^Niels Bohr Institute, Universitetsparken 5, DK-2100 Copenhagen, Denmark 
^OPMC, University of Geneva, Quai Ernest Ansermet 24, 1211 Geneva, Switzerland 



1 



The system is placed on a ring of circumference L and periodic boundary conditions 
are imposed: 

V't(o) = ^t(i), V(o)=V'(i). (1.3) 

The particle number operator is 

N = ! dxi^\x)i:{x), (1.4) 
Jo 

and the momentum operator P is 

p=l/ dx \^ d:,i)\x) i:{x) - il,\x)da,il>{x)Y (1.5) 
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These operators are integrals of motion: 

[H,N] = [H,P]=Q. (1.6) 

The interaction constant c > has the dimension of inverse length. The dimensionless 
coupling strength 7 is given by 

7=^, (1.7) 

where D is the particle density, 

= f (1.8) 

The model has been studied extensively. Its exact eigenfunctions, spectrum, 
and thermodynamics were obtained during 1960s, while the calculation of its correlation 
functions still remains a challenge (for a review of the field see, for example, Ref. |^]). 
We obtain in the present paper an exact expression for the local correlation function 

53(7) = {i^Hxfi^ixf), (1.9) 

where (• • • ) denotes the expectation value in the ground state |gs) of the system-^, 

(•••) = (gs|---|gs). (1.10) 

Let us formulate the main result: in the limit of the infinite number of particles taken 
at a finite density, 

D — > const 7^ 0, 00 as N,L ^ 00, (l-H) 
the correlation function (^]^) is expressed as 



53(7) 3 , , 7\ / 0^2 3e2e2 , 9e 



, e4-^+(l + 7^)4-2--^ + ^, (1.12) 

n3 27 72 V 2/ 7 7 72 ^ ^ 

where em are moments of the quasi- momentum distribution function (T{k) 

J dzz^aiz), (1.13) 



^Since the ground state is translationally invariant, 53 do not depend on x. The matrix elements, in 
particular, the expectation values, of local operators, like ijj'^ {xY'ip{xY , are often called form-factors. 
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and is the derivative of with respect to 7. The function cr{k) is the solution to 
the linear integral equation 

1 2aa(y) 1 

a{z) -— dy , , = — , (1.14) 



2-K y_i + {z- yY 2ti 

where a is an implicit function of 7 : 

a = 'J J dza{z). (1-15) 

The paper is organized as follows: In section |2| we consider the g-boson lattice 
model. This is an integrable lattice model whose continuum limit is the Lieb-Liniger 
model. The studies of the Noether currents in the g-boson lattice and the Lieb-Liniger 
model are performed in section |3|. In section ^ the g-boson lattice model is bosonized. 
In section ^ we write some contour-independent integral identities relating short- and 
long-distance properties of correlation functions in the g-boson lattice model. At long 
distances the results of section |^ are applicable. We thus get some non-trivial informa- 
tion about the short-distance properties of correlation functions. In section |^ we take 
the continuum limit of expressions obtained in section ^. This leads us to Eq. ( |1.12|) 
for 53(7). This equation and the contour- independent integral identities of section |5| 
are the main results of our paper. 



2 g-boson lattice model 



In this section a lattice regularization of the Lieb-Liniger model ( |1 . ![ ) is described. Lat- 
tice regularization is a natural tool to circumvent the short-range singularities discussed 



in section 3.2 and can, in principle, be done in many different ways. An integrable lat- 
tice regularization will be needed for our purposes. Several different regularization 
schemes preserving integrability are discussed in Ref. Q. We shall use the so-called 
g-boson lattice model [0, |8| as an integrable lattice regularization of the Lieb-Liniger 
model. This choice is motivated by its simplicity. 



In section p.l| the g-boson algebra is constructed for a given lattice site. In section 
|2.2| we present the Hamiltonian of the g-boson lattice model defined on a lattice with 
an arbitrary number of sites, M. In the continuum limit, defined by Eq. ( |2.19| ), this 
Hamiltonian becomes the Hamiltonian of the Lieb-Liniger model, Eq. ( |1.1| ) . In section 



2.3 we discuss the integrability of the g-boson lattice model. The generating functional 
of the integrals of motion is constructed using the formalism of the Quantum Inverse 
Scattering Method. By expanding this functional the explicit expressions for several 



integrals of motion are obtained. In section 2A we study the eigenfunctions and the 
eigenvalues of the integrals of motion in the g-boson lattice model. Finally, we apply 
the Hellmann-Feynman formula to one of the integrals of motion in section |2.5| with 
the aim of using the resulting identity in subsequent sections. 

2.1 g'-boson algebra 

Consider the operators -B„, b\ and Nn = nX satisfying the so-called g-boson algebra 

BnBi-q-^BlBn = l, (2.1) 
[Nn, Bn] = -Bn, [N^, bI] = bI (2.2) 
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The index n labels the sites of a lattice; it will be held fixed in this section, where we 
work with the operators for a given lattice site. The parameter g is a c-number and it 
is enough for our purposes to work with real q > 1. 

Define now a Fock space for the g-boson algebra given by Eqs. ( |2.1D and ( |2.2D . Let 
the basis states 1^-^)^, nij = 0, 1, 2, . . . of the Fock space be the states of the harmonic 
oscillator 

bn\mn)n = ml/'^lmn - l)n, ^Jj |m„)„ = (m„ + 1)^/^ |m„ + 1)„ , (2.3) 

where bn and 6„. ai^e the canonically commuting creation and annihilation operators 

[bn,bl] = l. (2.4) 

In this Fock space the operator A^^^ entering Eqs. ( |2.2|) acts in a manner identical with 
that of the operator 

Nn = bibn, (2.5) 

for which we have 



Nn\0)n = 0, Nn\mn)n = inn\mn)n, m„ = l,2,3, (2.6) 

The operators Bn and Bn entering Eqs. ( p.l| ) and (2.2) act in the Fock space as follows 

Bn\mn)n = [mn]y^\mn - l)n, Bl\mn) n = [mn + Ij^^ \mn + I) n , (2.7) 



where 



-2x 



7-2 ■ 



(2i 



It can be easily shown that the commutation relations (|2.1| ) and ( |2.2| ) are consistent 
with Eqs. (U) and Note that 



as 



1, 



(2.9) 



and therefore 

Bl bl and Bn b^ as 1. (2.10) 

It should be emphasized that the present paper deals with a representation, Eqs. (2.6) 
and ( |2.7D , of the g-boson algebra, but not with the algebra itself. All statements con- 
cerning the g-boson algebra should be understood as statements about this particular 
representation. One can see from Eq. ( p.5| ) that this choice of representation makes it 
possible to relate the operator Nn entering Eqs. ( |2.2D with the canonical boson oper- 



ators b\i and 6„. It is also possible to express the operators Bn and Bn in terms of b\i 
and bn : 



Br, 



'[Nn + l]q 



bn, 



"V Nn + l 



V Nn + l 

and give an alternative form of the commutation relation (2.1): 
The relations ( 2.11| ) and ( 2.12| ) can be easily proven using Eqs. (| 



(2.11) 



(2.12) 



and 
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We see from Eqs. ( p.5D and (2.11) that the relations between the g-boson operators 
Bn and B^ and canonical boson operators bn and 6„ are nontrivial. Furthermore, it 
follows from Eqs. (|2.7| ) that 



BlBn\mn)n = [^IT-nlg (2.13) 

and therefore that 

BiBn = [Nn]g = • (2.14) 

Thus, BnBn 7^ Nn, and is non-polynomial in terms of Bn and Bn- 

2.2 Hamiltonian of the g-boson lattice model 

In order to pass from the quantum mechanical (one-site) model considered in the 
previous section to a quantum field theory model we introduce a lattice with M sites; 
the index n in Bn, b\ and Nn labels the lattice sites. We impose an "ultralocality" 
condition on the operators Bn^, Bn^ and Nng by requiring that they commute at 
different lattice sites. The basis states of the whole lattice are constructed as the 
tensor product of the local basis states: 

|0) = ®ji^\0)j, \m) = ^ji^\mj)j. (2.15) 

The Hamiltonian for the (;-boson lattice model is defined as follows: 

1 ^ 

H, = -^Y.^BiBn+i + Bl_^,Bn-2Nn) (2.16) 

n=l 

with the periodic boundary conditions M+1 = 1 imposed. Here 5 is the lattice spacing. 
The factor d~'^ is introduced to ensure the proper continuum limit. One can check using 
the g-boson algebra, Eqs. ( ^.l| ) and ( |2.2D , that Hq commutes with the number operator 

M 



N = J2Nn- (2.17) 



n=l 



Since N is non-polynomial in terms of Bn and Bn, which can be seen from Eq. (|1|), 
Hq is non-polynomial in terms of these fields. It is non-polynomial in terms of the 
canonical lattice bosons bn and 6„ as well, now due to the non-polynomiality of the 
hopping term. Thus the model is interacting and the interaction is encoded in the 
deformation parameter q : if one takes g — > 1 the result will be a quadratic boson 
Hamiltonian on a lattice, describing nearest-neighbor hopping: 

1 ^' 

Hq^ —^Y.{blbn+i + bi^^bn-2Nn) as q ^ 1. (2.18) 

n=l 

The limit g ^ 1 is somehow trivial. It is much more fruitful to consider the following 
limit: let 5 — > 0, M — > oo, and q ^ 1, while L and c are kept constant: 

L = M5, c/2 = Hi6~~\ as <5 ^ 0, M ^ oo, k ^ 0, (2.19) 
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where k is related to q as follows 

q = e''. (2.20) 

The limit ( ^.19 ) will be called the continuum limit of the g-boson lattice model. The 
sums in the continuum limit are converted into integrals in the usual way 

M ,.L 

(5 V ^ / dx. (2.21) 

n=l 

For any g-deformed quantity [x\q, Eq. ( |2.8| ), the following expansion is valid 

^2 

[x\q = x-Kx{x-l) + —x{x-l){2x-l) + ..., K ^ 0, (2.22) 

3 

where q is related to k by Eq. ( p. 20 ). Therefore 



y^^^ = l-fA^n + ^iVn(5iV„+4) + ..., .^0. (2.23) 

We define the continuum boson fields ^{x) and tp^x) by 

= 8~^'Hl V(x) = (2.24) 

where 

x = 5n. (2.25) 
The fields ip\x) and il>{x) satisfy canonical commutation relations (|1.2|). One can easily 



check that the g-boson Hamiltonian (|2.16| ) becomes the Hamiltonian of the Lieb-Liniger 
model (1.1) in the continuum limit ( p. 19 ). 

Finally, we discuss the number of particles and momentum operators in the g-boson 
lattice model and their continuum limit. We have defined the number operator in the 
g-boson lattice model by Eq. (|2.17| ). The continuum limit ( ^.19|) of this operator is 
given by Eq. (|1.4| ). The momentum operator is defined in the g-boson lattice model as 
follows: 

^1 = -h J2(^nBn+l - i?l+li?n). (2-26) 
n=l 

Its continuum limit is given by Eq. (|1.5| ). To avoid any confusion, we note that Pq is 
not a generator of lattice translations, that is 

e-^^''fi„e^^« (2.27) 

However, its continuum limit (|l]^) is such a generator for the continuum model: 

e-'''^tp{x)e'''^ = ^{x + a). (2.28) 

It is due to the property ( |2.28| ) that we call Pq the momentum operator. We shall not 
need to know an explicit form of the true momentum operator P'^^'^ which generates 
the lattice translations in the g-boson lattice model. 
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2.3 Integrals of motion in the g-boson lattice model 



An infinite set of integrals of motion in the g-boson lattice model can be constructed 
using the Quantum Inverse Scattering Method. A detailed description of this method 
is given, for example, in Ref. and it was applied to the g-boson lattice model, in 
particular, in Refs. [0, |8|. We give in the present section a schematic description of the 
method, referring the reader to the Refs. [0, ^, ^ for details. The main purpose of the 
section is to obtain Eqs. (|2.44 )--( ^^ . 

The so-called L-operator for the g-boson lattice model is defined by 



LniX) 



xBi 



XBn 



where 



X 



-2k 



(2.29) 



(2.30) 



The L-operator (|2.29|) is a 2 x 2 matrix with the entries being quantum operators 



acting in the infinite-dimensional Fock space defined by Eq. (2.3). The Quantum 
Inverse Scattering Method is based on the existence of the intertwining relation for the 
L-operator: 

R{X, ^) (L„(A) Ln{^l)) = (L„(/i) L„(A)) R{X, fi) 



(2.31) 



with the i?-matrix defined by 

/ /(/^,A) 

i?(A,M) = 



V 











g-i 






q 





where 



/(A,/i) 



sinh(A ~ fx + k) 



gix,n) 



sinh(K 



sinh(A — /x) ' 'f*/ sinh(A — ^) 

Recall that the tensor product for 2x2 matrices a and b is defined as follows: 



(2.32) 



(2.33) 



an ai2 

021 022 



hi bi2 

&21 &22 



/ aiiftii 011&12 auhi a^bu \ 

011621 011^22 012^21 012^22 

021^11 021^12 022^11 022^12 

\ 021621 021622 022621 022622 / 



. (2.34) 



The monodromy matrix T(A) is defined as a matrix product of the L-operators taken 
over all lattice sites 



T(A) 



^(A) B{X) 
C{\) D{X) 



Lm{\)Lm-i{X)---Li{\). 



(2.35) 



The entries A{\), . . . ,D{X) of the monodromy matrix are quantum operators acting 
in the tensor product of the local Fock spaces over all sites of the lattice. Due to 
the relation (2.31) and the commutativity of the entries of the L-operator ( |2.29| ) at 
different lattice sites one has the intertwining relation for the monodromy matrix 



R{x, fj) (r(A) ® r(^)) = (r(^) ^ r(A)) r{x, ^). 



(2.36) 
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Equation ( p.36| ) defines 16 commutation relations for the operators entering the mon- 
odromy matrix. We write explicitly those relations which we shall use in deriving the 
Bethe equations in section 



qAiX)Bi^l) = /(A, fi)B{i^)A{X) + g{^l, X)B{X)A{fi), (2.37) 
qD{X)B{n) = /(^, X)B{fi)D{X) + g{X, /i)S(A)I)(/x), (2.38) 
[B{X),Bif,)]=0. (2.39) 

The transfer matrix r(A) is defined as the trace over the matrix space of the mon- 
odromy matrix 

r(A) = Trr(A) = A{X) + D{X). (2.40) 
It can be proven (see chapter VI of Ref. Q for details) that for any A and 

[r(A),r(/.)]=0 (2.41) 

which implies that r is a generating function of the integrals-of-motion of the problem: 
expanding t(A) in A one gets a set of commuting integrals-of-motion Im- This set can 
be chosen in many different ways since any analytic function of t(A) can play the role of 
the generating functional. We, however, impose an additional very restrictive locality 
condition on these integrals by requiring them to be written in the following form 

M 



SY.J^^'Hn), (2.42) 



ra=l 



where the operators Jt^^ {n) act nontrivially in m neighboring lattice sites only. The 
subscript r appears in jT^\n) since these local operators are r-components of the 
corresponding Noether currents, which will be discussed in section |3|. To get the set 
{Im} we introduce the variable 

C = (2.43) 



and consider the expression 

1 



(2m)! dC^" 



In [C^'r(C)] 



m = 1,2,3,.... (2.44) 



We assume that A is real, therefore C is real and nonnegative. The local operators 
J^^\n), jP{n) and J^^\n) are 

jy\n) = ^X^BlBr,+i (2.45) 

jj^\n) = (l - y) {BiBn+2 - ^^BiBiBn+A+, - BlBl^^Bn+A^^ 

(2.46) 
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and 



J^''\n) = {l-X^ + y) {BiBn+^ - x^BiBiBn+iBn+2 - x^BiBl^^Bn+iBn+z 



■X 



BiB^BlBn+lBn+lBn+l 



3-3x2 + x^^ 

+ X'^BlBlBl_^_^Bn+lBn+lBn+2 + X'^^i^i+lfiI+l^n+l-B„+2-B„+2 



+ x^BiBi^^Bi^2Bn+iBn+2Bn+3^ (2.47) 



respectively. To calculate 53(7), Eq. (|1.9|) , we will need J^r and ; the explicit 

(3) (m) 

expression for jTr is displayed in order to illustrate how the complexity of jTr grows 
with increasing m. 

The integrals Im are non-Hermitian, they contain both real and imaginary part. 
Using the involution 

[r(C)]^ = riC') 



K2). 



it may be shown that 

[/^,7„]=0 for any m,n. 
For the integrals of motion Im, we use the following notation 

I~m ~ Im' Ttl = 1, 2, 3, ... , 

and for the local operators Jr 

jy'^Hn)^[J^^Hn)]\ m= 1,2,3,. 



Using the involution ( 2.4^ ) one gets for 

1 (i^™ 



(2m)! dC^™ 



In [e^TiQ-^)] 



m = 1, 2, 3, , 



(2.48) 
(2.49) 

(2.50) 
(2.51) 

(2.52) 



A set of common eigenfunctions of J^, rn = ±1, ±2, ... is constructed in section 2.4. 

The local operators Jt^^ {n) generated by Eqs. ( 2.42| ) and (2.44) are polynomial in 
Bn and Bn, while the one-site number operator A^^ is non-polynomial in these variables, 
as was noticed below Eq. ( p.l4| ). Therefore, the number operator A^, Eq. ( p.l7| ), cannot 
be expressed as a finite linear combination of the integrals of motion Im, Eq. ( |2.44| ). It 
is, however, clear from the structure of Jr'^^n) that 



m = ±1,±2, 



(2.53) 



Indeed, A^ commutes with any monomial containing an equal number of the creation 
and annihilation operators, B^ and B. It follows from Eq. (2.53) that 

[N,lm]=0, m = ±l,±2,.... (2.54) 

To stress that A^ is one of the integrals of motion, we shall use the following notation 



A^ = /n 
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1 



(2.55) 



The Hamiltonian ( p. 16 ) can thus be written as 



^1 = ^(^1 + - 2x'^o), (2.56) 



and the momentum operator ( |2.26| ) as 

i 1 



Pi = -i^ — ih-I~i)- (2.57) 



2.4 Eigenfunctions and eigenvalues of the integrals of mo- 
tion in the g-boson lattice model 

We have constructed in section a set of integrals of motion Im , Eq. ( |2.44| ) , of the 
g-boson lattice model on a lattice with an arbitrary number of sites, M. We find 
in the present section their common eigenfunctions and their eigenvalues using the 
Algebraic Bethe Ansatz technique an important ingredient of the Quantum Inverse 
Scattering Method. The eigenvalues of all Im are defined in the A^-particle sector by 
N parameters called quasi-momenta. These quasi-momenta are the solutions of the 
system of nonlinear equations called the Bethe equations. Following the presentation 
of Ref. 1^, Chap. I, we find from the analysis of the Bethe equations the ground 
state quasi-momentum distribution and study its properties in the limit defined by 
Eq. 1 ^ ). 



The cornerstone of the Algebraic Bethe Ansatz method is the fact that the vacuum 
|0) defined by Eq. ( |2.15| ) annihilates the C(A) entry of the monodromy matrix ( 2.35| ) 
and is an eigenfunction for the A(X) and D{X) entries: 

A{X)\0) = e^'^^\0), D{X)\0) = e-^'^^\0). (2.58) 

Since is a good quantum number, one can work in the A-particle sector. Define a 
set of states by the following formula 

N 

|V^^(Ai,...,A^)) = n^(Ai)|0). (2.59) 

i=i 

These states (often called the Bethe states) are eigenfunctions of the transfer matrix 
(2.40) and, hence, of the integrals of motion Im, m = ±1, ±2, . . . , if the parameters 



Ai, . . . , Am satisfy a system of coupled nonlinear equations (called the Bethe equations), 

^'""'^ = n TTTXT' J = h...,N. (2.60) 

fc=l •'^ 

The eigenvalues 9^ of the transfer matrix acting on the Bethe states ( |2.59| ) are given 

by 

t(A)|V'7v) =^7v(A,{A,})|7/;7v), (2.61) 

N N 
q^ONiX, {A,}) = e*^^ n /(^' ^^■) + n f(^J^ (2.62) 
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The eigenvalues of the integrals of motion /„ 



m 



±1,±2, 



acting on the Bethe 



states ( |2.59| ) can be obtained by acting with the representations ( 2.44| ) and ( |2.52| ) onto 
I'i/'Tv) and using Eqs. (|2l^ ), (|]6l|) and (|^. The calculations are tedious, while the 
final result is surprisingly simple: 



-2\m\ 



1 ^ 



m 



±1,±2, 



(2.63) 



It is useful to mention that, in order to calculate the correlation function (L9) we shall 
need to know the spectrum of the integrals I±i and I±2 only. For \m\ > 2 we did not 
obtain Eq. ( 2.63| ) analytically. Instead, we simply checked that it is correct for some 
given values of A^, M and m using the Mathematica package. 

It will be convenient to use instead of Ai, . . . , \n a set of quasi-momenta pi. 



Pj 



J 



~ 2 

Written in these variables, the Bethe equations 

N 



n 

k=i 



Pk) + ^K\ 



k=i ^^"[^ 



{Pj -Pk)-ii^] 



are 



J 



• ,PN ■■ 

(2.64) 
(2.65) 



Using Eqs. ( p. 56] ), ( |2.63| ) and (2.64), one gets for the eigenvalues of the Hamiltonian 



N 



2P± 

2 



(2.66) 



We now discuss some properties of the Bethe equations necessary to identify the 
ground state of the model and to take the limit ( 2.80| ). The analysis will be very 
similar to that one carried out for the Lieb-Liniger model in Ref. [^j. Chap. I, so we 
omit several long proofs, referring the reader to Ref. |9|1 for details. 



(i) All the solutions pj of the Bethe equations ( 2.65| ) are real. The proof is the same 
as that given for the Lieb-Liniger model in Ref. [^, page 11. 

(ii) It follows from Eq. (|]6|) that En is a periodic function of the quasi-momenta 
Pj with period 2tt. We shall work with quasi-momenta lying in the interval 



-vr < Pj < vr. 



iV. 



(2.67) 



The condition ( |2.67| ) will be assumed in all subsequent formulas, 
(iii) We write the Bethe equations ( |2.65 ) in the logarithmic form 



N 



Mpj+Y,e{pj 



Pk) 



27r 



k=l 



Uj + 



A^- 1 



where the parameters rij take arbitrary integer values: 

rij is an arbitrary integer, J = 1, 
From now on we shall work with odd N : 

N is odd. 



(2.68) 



(2.69) 



(2.70) 



11 



The function 9{p) is 



0(p) = iln!!^4^^1±|4, 0{-p) = -e(j,). (2.71) 



sm (IK 



The derivative of 6{p) is positive, 

X ^ ^sin(2iK) ^ sinh(2/t) ^ ^ 

2 sin(iK; + ip) sin(iK - ^p) cosh(2K) — cosp 

therefore 9{p) grows monotonously in the interval — vr < p < tt (recall that the condition 
( 2.67D is assumed). 

(iv) For any set {nj}, Eq. ( |2.69| ), there exists a uniquely defined set {pj} of solutions 
of the Bethe equations ( p.68| ) . The proof is the same as that given for the Lieb-Liniger 
model in Ref. Q, page 12. 

(v) We write, using Eq. ( |2.68|) , 

N 

M{pj -Ps) + Yl ^^^Pi - - ^(P' - ^'fc)] = - ^s), j,s = l,...,N. (2.73) 

k=l 

Since 0'{p) > 0, Eq. (|2.72| ), the left hand side of the Eq. (p.73|) is a monotonically 
growing function of the parameter pj — p^. Therefore, if rij > then pj > ps] if 
rij = Us then pj = pg. We have thus shown that the set of quasi- momenta {pj} is 
uniquely characterized by the set {nj}, and vice versa. 

(vi) Like for the Lieb-Liniger model (Ref. |^], page 14) the energy functional ( 2.66| ) 
taken on the sets {pj } of the solutions of the Bethe equations has the minimum in the 
sector with the fixed number of particles, A^, if nj take the values 

nj = -N + j, j = l,...,N. (2.74) 



The Bethe equations (|2.68| ) for the ground state are, therefore, 

N 

Mp, + G{pj - Pk) = 27r [ j - 1 - ) , J = 1, . . . , AT. (2.75) 

k=l ^ ^ 

The ground state is obviously non-degenerate. 

(vii) It is obvious that the set {pj} of the solutions of Eq. ( p. 75 ) is symmetric with 
respect to zero. It follows then from Eqs. ( ^.631) and. (|2.6^ ]) that the eigenvalues of 

m = ±1, ±2, . . . corresponding to the ground state wave function |gs) are real, and 

{Ira) = {I~m), m = 1,2,3,... (2.76) 

where I-m is defined by Eq. ( ^.50]) . 

(viii) The ground state wave function is translationally invariant. This implies that 
the ground state average (j7T™^(n)) is n- independent. Using this property, one gets 
from Eq. (|2l^ ) 

{Im) = 5M{J^^^). (2.77) 
Comparing this with Eq. ( |2.76| ), one arrives at 

(J^™)) = (jM)^ (2.78) 
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where J'j ™^ is defined by Eq. ( p.51|) . We recall that the size of the system, L, is the 
product of the number of the sites, M, and the lattice spacing, 6, Eq. (2.19), therefore 
the particle density, Eq. ( |1.8| ), can be written as follows 

D = ^ = {j(0)). (2.79) 

We are interested in the ground-state properties of the g-boson lattice model in the 
limit 

^ oo, M ^ oo. (2.80) 

We introduce the quasi- momentum distribution function p{pj) by means of the follow- 
ing identity 

N N 

Y. = mY, P(.Pj)iPj+i - Pj), (2.81) 
j=i i=i 

where 

1 

M{pj+i -Pj)' 

The quasi- momenta pj fill the symmetric interval —A < pj < A and one has pj+i —pj ~ 
(Ref. Q, page 14). For an arbitrary function f{pj) one has 

N .A 

V/(p,)^m/ dpp{p)f{p). (2.83) 

The parameter A plays a role analogous to that of the Fermi momentum: all states 
with \p\ < A are occupied, and all states with |p| > A are empty. The value of A is 
defined by the normalization condition 



piPj) = T77Z —■ (2-82) 



N 
M 



A 

dpp{p). (2.84) 



The Bethe equations (2.75) become in the limit ( |2.80| ) a linear integral equation for 
the ground state quasi-momentum distribution p{p) 

1 1 

P{p) ~ ^ J ^ dpK{p - p)p{p) = — , (2.85) 

with the kernel K{p) given by Eq. ( 2.72 ), 

cosh(2K) — cosp 

Equation ( p.85| ) is often called the Lieb equation. 

Having defined p{p) and A one can easily get the ground-state expectation values 
of the integrals of motion, (/m), in the limit ( p. 80 ). In particular, 

(A) _ . 2 



M 

and 

M ^ 



X I dp p{p) cos p, (2.87) 

A 



(l-y) J ^dpp{p) cos 2p. (2.^ 
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The ground-state energy of the system can be obtained from Eq. ( 2.56| ) or (|2.66| ): 



En 



M 



^ J ^dpp{p) sm^|- 



(2.89) 



We shall make use of the properties of (Im) in section 6.1 



2.5 Hellmann-Feynman theorem 

In this section we use the Hellmann-Feynman theorem to derive an identity for the 
ground-state average of some local operator in the g-boson lattice model. This identity 
is given by Eqs. ( |2.92| ) an d (|2.96| ) below, and will be used in section |6.2| . 

It was shown in section 2.4 that the eigenfunctions ( p. 59 ) of the g-boson Hamiltonian 
(2.16) are the eigenfunctions for the integrals of motion Eq. ( 2.44| ) and their Hermitian 
conjugate Eq. ( ^.50|) . Since /i|gs) = /||gs), Eq. ( 2.76 ), one has 



(2.90) 



where (• • • ) denotes the ground state average and q is the deformation parameter 
defined by Eq. (^.1|). Equation ( 2.9C| ) is known under the name of the Hellmann- 
Feynman theorem. The explicit expression for Ii via local fields and is given by 
Eqs. (|2.42| ) and ( 2.45 ). Using the translational invariance of the ground state, implying 
that 



{Bl.-.BjJ 



(B 



t 

'ji+k ■ ■ 



where k is an arbitrary integer, equation ( ^.90| ) can be written as follows: 



Mi —{B\B. 
dq 



(2.91) 



(2.92) 



The dependence of the operators Bj and b] on the parameter q can be analyzed 
using the representation (2.11|). One gets from this representation 



Introducing the notation 
and combining Eqs. ( 2.93| ) and (2.94) one finds 



x + 1 
72^+2 - 1 



(2.93) 
(2.94) 

(2.95) 



d 



{B]Bj+i) = B]g{q,N,)Bj+i + B]g{q, Nj+^)Bj+^. 



(2.96) 
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3 Noether currents 



Noether's theorem is an important ingredient of Quantum Field Theory. In short, 
it states that symmetries imply conservation laws. More precisely, if the action of 
a system is invariant under an infinitesimal transformation of the fields, then there 
exists a function of these fields whose divergence is zero. This function is called the 
Noether current associated with the symmetry. One can say more about this function 
if the symmetry transformation leaves the Lagrangian (or, even better, the Lagrangian 
density) and not just the action invariant [p!o|]. 

For integrable models, however, it is often more natural to work within the Hamil- 
tonian formalism rather than within the Lagrangian one. This is due to the fact that 
explicit expressions for the integrals of motion are usually known in these models. 
In case of the g-boson lattice model the local integrals of motion are generated by 
Eq. ( |2.44| ) and explicit expressions for like (|2^)-(|2^, can be, in principle. 



written down up to an arbitrarily large m. The operator Ji^ can be recognized as 
the imaginary-time component of the conserved current J^^^ . By calculating the com- 
mutator of J'r™^ with the Hamiltonian one gets the space component c/i'"^ of the 
corresponding conserved current. The continuity equation for will be extensively 
used further derivations. 



In section 3.1 we recall briefly the notion of the Noether currents in the Hamiltonian 



formalism. Noether currents in the Lieb-Liniger model are considered in section 3.2 
We demonstrate the problem of making an unambiguous definition of higher integrals 
of motion in this model. It is because of this problem that we are working with the q- 
boson lattice regularization (described in section ^) of the Lieb-Liniger model. Noether 
currents in the q-boson lattice model are considered in section |3.3|. Finally, in section 



3.4| we subject the g-boson lattice model to a local gauge transformation and study the 



behavior of the Noether currents under this transformation. 

3.1 Noether currents in the Hamiltonian formaUsm 

To begin with, let us introduce some notation. It is often more convenient in Quantum 
Field Theory to work with the Euclidean (imaginary) time r rather than with the 
Minkowski time t : 

r = it. (3.1) 
We denote by 5^,^ the metric tensor for two-dimensional space-time theories^: 

_ J diag(l,— 1) Minkowski a o\ 

" \ diag(l, 1) Euclidean ^"^^^ 

The summation is performed over the contracted indices, 

a^foM = ^a^5/^ (3.3) 

and the rules for converting between covariant and contravariant indices are 

a;. = 5M.«', a^ = 9'"'a,, g^,g-- = 5Z. (3.4) 



^We choose Qtt = — 1 and gxx = 1 in Minkowski space-time. 
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The Hamiltonian H oi a theory can be defined as the generator of time-translations 
for an arbitrary operator O : 



O{x,t) = e''^O{x,0)e-''^. (3.5) 

We assume that the Hamiltonian is time-independent: 

dH dH , , 

— = — = 0. 3.6 

dt dt ^ ' 



Then the equation of motion for 0{x,t) following from Eq. (|3.5D is 

^OM{x,t) = j^OM{x,t) = i[H,OM{x,t)]. (3.7) 

We note that all the dependence on t of the operator OMix,t) is explicitly given by 
the evolution operator e**^, Eq. (3.5), and there is no need to distinguish the partial 
d/dt and full d/dt derivatives. The subscript M is used in Eq. (3.7) to indicate that 
it is written in the Minkowski time t. Written in the Euclidean time r, Eq. ( |3.1| ), the 
equation of motion (|3.7|) takes the form 



^Oe{x, t) = ^Oe{x, t) = [H, Oe{x, t)], (3.8) 

where Oe{x,t) = Om{x, t). We will mainly work with the Euclidean time and therefore 
we drop the subscript E to shorten notation. We shall also drop one or both of the 
arguments x,t in 0{x,t) when they can be recovered from the context. 

Impose a locality condition onto the operator 0{x) : suppose that this operator, 
taken at a point x, depends on the basis fields tp and ip^ (and their derivatives) taken 
at this point exclusively. Equations ( |3.1^^ and ( |3.16| ) provide us with examples of the 
operators of such type. If, in addition, dx 0{x) is an integral of motion, 

[H, [ dxO{x)] = 0, (3.9) 
Jo 

then 

[H,0{x)] = -dMx). (3.10) 

The symbol J'x is called the j;-component of the conserved (Noether) current. The 
operator O itself plays a role of the r-component of the Noether current: 

Jr{x) = Oix). (3.11) 

Indeed, combining Eqs. (^), ( p. 10 ), and ( |3.11[ ) one gets 

drJr + d^Jx = drJ^ + d^J"" = (Euclidcau) (3.12) 

which is the continuity equation for the conserved (Noether) current {Jt,<Jx)- We 
thus derived Noether's theorem within the Hamiltonian formalism: to every integral 
of motion, Eq. (|3.9| ), there corresponds a conserved current, Eq. ( 3.12| ). It should be 
stressed that the time derivative drJr is defined by the right hand side of Eq. ( |3.8| ); 
having drJr-, one can use Eq. ( 3.12 ) to calculate Jx- 
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3.2 Noether currents in the Lieb-Liniger model 



We now consider the Lieb-Liniger model ( |1.1| ). The local density J\f{x) of the number 
operator Eq. (|1.4|) is 

M{x) = J^^\x) = ip^x)ip{x). (3.13) 



(3.14) 



Commuting Af{x) with the Hamiltonian ( |l.lj ) one gets 

[H,M{x)] = -d,j('\x), 

where 

jP{x) = [d,4^\x)]i;{x)-i;\x)d^i;{x). 



(3.15) 



Another conserved current is the current associated with the local density V{x) of the 
momentum operator (|1.5P 



(3.16) 



V{x) = ^ I ^cl:^p\x) ^(x) - ip'^ {x)^x^p{x)^ . 

Commuting V{x) with the Hamiltonian ( |1 . 1| ) one gets 

[H, V{x)] = -d-^M{x) (ill-defined, see Eq. (p5|)), 

where 



(3.17) 



M{x) = ^^ dli;\x) ip{x)-2 d^ip^x) d^ipix) + ^p\x)^l^/J{x) - 2ctlj^ {xf^Pixfy 

(3.18) 

Let us discuss Eq. ( 3.18D in more detail. In getting this expression, one necessarily 
introduces objects like ^/^^(x)9^V(^) (this is clearly seen from Eq. ( 3.17] )) and we have 
assumed that such objects are well-defined. This is, however, an incorrect assump- 
tion. To show this we use results from Ref. [^l|. Consider the ground state average 
{Tp^{x)Tp{0)). When x ^ one can expand ip^{x) in the Taylor series: 



{i^HxMO)) = (V'^(O)V(O) + ^9,^t(e)^(o) 



1! 



. (3.19) 



e=0 



All the terms written explicitly on the right hand side of Eq. ( p. 19 ) can be found in 
Ref. lH: 

{^\0)m)=D, (3.20) 



where D is the average density, Eq. (pT^), 

{^,^|;He)m) 

and 



const (7) D^, 



where const (7) depends on 7, Eq. ( |1.7| ), exclusively. Finally, 



(5fV^(e)^(0)) 



sgn(e) X const(7) 



(3.21) 
(3.22) 

(3.23) 
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where 

sgn(e) = \ _^ g ^ Q • (3.24) 



Equation ( 3.23| ) is of crucial importance. It shows that 



[dli;^x)]'ilj{x + e) / [d^ip^ {x)]ip{x - e) when e ^ (3.25) 

so one should indicate explicitly how the point-splitting procedure is performed when- 
ever writing [d^ilj^ {x)]iIj{x). The same is true for other operator products containing 
d^, (in particular, Eqs. ( p.l7|) and ( p.l8|) need such a prescription) and more generally, 
for operator products containing the derivatives 9" with n > 3. We perform the point- 
splitting procedure by putting a system on a lattice and we discuss the corresponding 



Noether currents is section 3.3 



3.3 Noether currents in the g-boson lattice model 



It follows from the results of sections and 3.1 that the g-boson lattice model contains 
an infinite hierarchy of Noether currents. Since this model is a lattice model, all 
currents are well-defined. This is a major advantage as compared to the Lieb-Liniger 
model, which, as it was shown below Eq. (3.18), suffers from short-range singularities. 
We obtain in the present section various relations between conserved currents in the 
g-boson model. These relations will be exploited in sections ^ and |5|. 

The lattice version of the continuity equation ( |3.12 ) for a conserved current = 

{Jr,Jx) is 

d 1 

■7rMn,r) + - [Mn,T) - Mn - 1,t)] =0 (3.26) 



dr 

where, according to Eqs. 
defined as follows: 



and ( ^.11[) , the derivative of jTr with respect to r is 
d 



dr 



(3.27) 



The Hamiltonian ( |2.16 ) commutes with the total number of particles, Eq. ( [2.171 ). The 



r-component, (n\ of the corresponding local current is given by Eq. ( |2.55| ): 

J,(0)(n) = iiV„. 





(3.28) 



We calculate the commutator [i/^, A'n] with the help of Eq. (2.2), then substitute the 
resulting expression into Eq. ( 3.26| ) and get 



Using Eq. ( |2.45| ) one can write Eq. ( 3.2S| ) as follows: 

jf\n) = -^^{j^\n)-J^-^\n) 

In the continuum limit ( |2.19| ) equations ( |3.28| ) and ( |3.29| ) become 

jf^{x) =^\x)^{x), 

jf\x) = kV'^(x)l ^(x) - V^(x)5,V(x), 



(3.29) 
(3.30) 

(3.31) 
(3.32) 
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thus reproducing the expressions ( |3.13| ) and (|3.15| ) for the r- and x-components of the 
conserved current j'^^^ in the Lieb-Liniger model. 

The a;-component of the current J^^\n) can be calculated in very much the same 
way as ji^\n). Indeed, the r-component of J^^^n) is given byEq. (HH). Substituting 
this expression into Eq. ( 3.26| ) and calculating the commutator [Hq, Jr^\n)\ with the 
help of Eqs. Q, Q, and ( p^ ), one gets 



jP{n) 



52 L 



+ BiBn 



(3.33) 



where x is defined by Eq. ( |2.3[1| ). For the x-component of the current ^•'(n) one 
gets, taking into account the relation ( |2.51| ), 



(3.34) 



Having Eqs. (|!3|) and { ^I^ one can easily calculate the operator ^A, Eq. ( |3.17| ), 
in the g-boson lattice model. The operator A4 is defined unambiguously on the lattice 
by equation ( 3.26| ) : 



1 



[M{n) - M{n - I)], 



where 'P{n) is the density of the momentum operator ( |2.26| ): 



n) 



(3.35) 



(3.36) 



(recall that the local densities on a lattice are defined according to Eq. ( 2.42D ). The 
resulting expression for ^A(n) is 



Min) 



i 1 
'26^ 



(3.37) 



In the continuum limit ( |2.19| ) this expression transforms to Eq. ( |3.18| ). 
Finally, we note the following property of ^ji™^ (n) 

[N,J^"''\n)]=0, m = 0,±l,±2,.... (3.38) 

The proof of this expression is the same as that of Eq. ( ^.53|) . Equations ( 2.53| ) and 
(3.38) play an important role in the bosonization procedure for the g-boson lattice 
model, carried out in section i.4. 



3.4 Local gauge transformations in the g'-boson lattice 
model 

The operator in the g-boson model generates the global U{1) rotation (it is often 
called the global U{1) gauge transformation) of the fields Bn and Bt- 

e-^'t'^ Bne'^^ = Bne'^, e''^^ B^e''*'^ = B^e''^ . (3.39) 
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The T-component of the Noether current (n) is invariant under this rotation: 

g-i0iV^M(^)e*0iV = ^M(„)^ m = 0,±l,±2,.... (3.40) 

Whenever one works with a system possessing a global symmetry transformation, it 
is useful to extend this transformation to a local one. The local gauge transformation 
coming from the extension of the global U{1) symmetry is performed by a unitary 
operator 



Qe = exp 



M 



n=l 



(3.41) 



where J'T^\n) is the density operator ( 3.28 ), the parameter e is an arbitrary real 
number, and r]{n) is an arbitrary function of the lattice coordinate n. The ^-boson 
fields Bn and Bn (called "matter fields" in Quantum Field Theory) are transformed 
by the operator ( |3.41D as follows 



(3.42) 



For r][n) = const, Eq. ( |3.42D reproduces the transformation law (3.39). 

The r-component of\7^™)(n) is not invariant under the action of Q^; its evolution 
is described by the operator 

J^"'\n,e) = Q;'j^'^\n)Q,, m = 0, ±1, ±2, . . . . (3.43) 

For the case when the function r/(n) is a linear function of n, 

r](n) = a5n + h (3.44) 

one can get the following representation of the right hand side of Eq. ( |3.43| ) : 

= e^^'^''™J,(")(n), m = 0,±l,±2,.... (3.45) 

For and jj^^ this expression can be checked by applying to Eqs. ( |2.45| ) , 

( 2.46 ), and ( 2.47| ), respectively. By investigating the structure of Eq. ( ^.44 ) one gen- 
eralizes this calculation to the case of an arbitrary m. It follows from Eq. ( 3.45| ) that 



ia6mjj"'\n). 



e=0 



We shall need the ground state expectation value of this equation 



(3.46) 



(3.47) 



e=0 



We will also need to know the action of onto the x-component of the current 
J^^^\n), Eq. ( 13.291) . Assuming that r]{n) is the linear function of n, Eq. ( 3.44 ), one 
gets from Eqs. ( ^ ) and ( ^IB]) 



e=0 



(3.48) 



20 



Taking the ground state average and using Eqs. (2.77) and (|2.7S|) one arrives at 



e=0 



— I7W\ = 2m (Ji) 



We define now the currents j*^™"^ by the formula 



m = 0,±l,±2,..., 



(3.49) 



(3.50) 



The continuity equation for j^™) is 

a^^jM^o. (3.51) 
The e-dependent current ji™''(n,e) is defined in the same manner as ji^\n,e) : 



m = 0,±l,±2, .... 



It is obvious that 



e=0 



e=0 



One can, therefore, write Eq. ( 3.46| ) as follows 



iadmjj^'^n), m = 0, ±1, ±2, . . . . 



e=0 



In particular. 



0. 



(3.52) 
(3.53) 

(3.54) 
(3.55) 



e=0 



4 Bosonization of the g^-boson lattice model. 

In this section we discuss an effective field theory describing the low-energy properties 
of the g-boson lattice model. This effective field theory is the free boson theory in 
one space and one time dimension, studied in great detail in many review articles and 
textbooks, for example in Refs. [|l2|, O]. Some of its basic properties are reviewed briefly 



within the coordinate space formulation of section 4.1, the others are given within the 



momentum state formulation of section 4.2. All the correlation functions of the free 



boson theory can be calculated explicitly, making it possible to classify the operators 
of the theory according to their anomalous dimensions. This is the subject of section 



In section [4.4| we perform the so-called bosonization procedure: we establish the 
correspondence between the operators of the microscopic theory (the q'-boson lattice 
model) in the low energy limit, and the operators of the free boson theory (shortly, we 
take the bosonized limit of the microscopic operators). We continue the bosonization 



procedure in section 4.5 where we express the c-number coefficients in the bosonized 



representation of microscopic operators via {Im), studied in section |2|. 
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4.1 Free boson theory in coordinate space 

In the present section we review briefly the free boson theory in the coordinate space 
formulation. The Lagrangian density in Minkovski space is 



1 



M 



2ttK 



1 



and the corresponding Lagrangian density in the Euchdean space is 



1 



27rK 



-{dr(i)f + v{d. 



(4.1) 



(4.2) 



(the rules for converting between Minkovski and Euclidean spaces are given in section 
|3.l| ). The parameter v has the dimension of velocity, the parameter K is dimensionless, 
it is often called Luttinger parameter. The action Sm in Minkovski space is defined as 



Sm 



dt dx C 



M- 



(4.3) 



The model is placed on a ring of circumference L, the boundary conditions are discussed 
in the paragraph below Eq. ( f4.ll ). To define the Euclidean action Se note that the 
weight function e**^*^ , which appears in the functional integral formulation of the theory, 
oscillates in Minkovski space, while in Euclidean space it should decay rapidly and is 
conventionally written as e~^^. Therefore 



and 



Se 



2'kK 



-Se 



dr dx 



(4.4) 



(4.5) 



The field is a free massless real scalar boson field. An equation of motion for this 
field is the wave equation: 

d'^(l) + v'^dl(l) = 0. (4.6) 
The canonical momentum 11 (the field canonically conjugated to (f)) is 

dCM 1 



n 



-dt 



-dr(, 



ddt4> ttKv ttKv 
It is often convenient to work with the field 9{x), defined as follows: 

u{x) = -d.eix). 

vr 



(4.7) 



(4.8) 



We shall switch freely between dr4', dxO, and 11 in subsequent formulas. Being quan- 
tized, the fields (p and 11 obey canonical equal-time commutation relations 



[n(x),(/)(y)] = -i6{x - y) 
The Hamiltonian H of the system has the form 

H = I dxHix), nix] 



^:K[nnix)f + ^[d.cPix)f 



(4.9) 



(4.10) 



22 



where the symbol :: stands for normal ordering, discussed in section |4.2| . Another 
important ordering prescription is the time-ordering T. For any two boson operators 
A{t) and B{t') it is defined as follows 



TA{t)B{t') 



A{t)B{t'), t>t' 
B{t')A{t), t<t' ■ 



(4.11) 



We assume that T acts on all the operators standing on the right. 

Boundary conditions are an important issue in the theory. We impose periodic 
boundary conditions on the operators dx4> and dxO : 



(4.12) 



Since the Hamiltonian (4.10) contains dx(p and dxO, rather than these fields themselves, 
one has 

(/.(x + L) = 0(x) + ai, B{x + L) = e{x) + 02, (4.13) 
where ai and 02 are some operators. Their properties are discussed in section 



4.2 Free boson theory in momentum space 

In this section we give the momentum space representation of the free boson theory 



(4.10). This representation is a convenient starting point for calculating the correlation 
functions of the theory, like those considered in section 4.3 . 
The momentum space representation of the (j) and 9 fields is 



(p[X, t) = (pQ + TTVtK — — TTX 

Ij Ij 



^ E ^ " sgn('Z)e-''^- (6te™*l'?l + 6_,e-*l'^l) (4.14) 



1/2 



and 



^x,t) = eo-TTvt +7rx- + - 2^ 



KL 



2lT 



qLK 



1/2 



e-^«^ (ftje™*!^! - 6_ge-™*l«l) , (4.15) 



respectively. The momentum space representation of the Hamiltonian (4.10) is: 

1 



H 



TTV 

2L 



K 



{N - NqY + KJ' 



(4.16) 



The summation index q in Eqs. ( 4.14 )- (|4.16| ) runs through the following set of values: 

27r 



—J, j = ±l,±2,.... 



(4.17) 



The operators bq (bq) are boson creation (annihilation) operators obeying canonical 
commutation relations 

[bq,bl,]=6qq,. (4.18) 
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Thus, the last term in the right hand side of Eq. ( [4.16| ) represents a set of decoupled 
harmonic oscillators with the frequencies v\q\. The first three terms on the right hand 
side of Eqs. (|4l4D and ( p^ ) give the so-called "zero-mode contribution". They all 



commute with bq and bq] the only nontrivial commutation relations between themselves 
are 

[J,(l)o] = -i, [N,9o] = i. (4.19) 



Comparing Eq. ( 4.13| ) with Eqs. ( 4.141) and (|4.15D , one can easily express ai and Q2 



via J and — Nq (No is a c- number). The normal ordering symbol :: standing in 
Eqs. ( 4.10| ) and ( [4.26| ) means that in any given monomial one should place the creation 



operators bq to the left of the annihilation operators bq. 

We denote the ground state of the theory as |0)); the ground state expectation value 
((O)) of an arbitrary operator O is ((OlOlO)). We use the symbol ((• • •)) for the ground 
state expectation value of the free boson theory in order to distinguish it from the 
ground state expectation value (• • • ) of the microscopic model, Eq. (|l.9|) . One has 



bk\0)) = for A; = ±l,±2, ... (4.20) 

and 

JI0)) = 0, iV|0)) = iVo|0)). (4.21) 

Finally, we define the action of the operators 0o and 9q on |0)). We do this taking the 
exponent of these operators^: the states e*™''^''|0)) and e*"^"|0)) are non-vanishing states 
orthogonal to the vacuum state: 

(^/^e'rn^o^) = 5mo, ((e*"'° )) = .5„o, (4.22) 
where 5mim2 is the Kronecker 6 symbol 

j I mi =1712 I. 
[0 mi ^ 7712 

Any operator within the free boson theory we shall work with, will contain the fields 
00 and ^0 in the form e*™"^" and e*"^" exclusively. 

4.3 Correlation functions and spectrum on the anomalous 
dimensions in the free boson theory 

The correlation functions of the free boson theory, Eq. (|4.10| ), are known explicitly 



and are given in many textbooks on one-dimensional physics p^ , 13 1. They can be 
calculated, for example, starting from the momentum space representation, discussed in 
section The knowledge of the correlation functions provides us with an efficient tool 
for classifying the operator content of the theory. We classify the operators according 
to their anomalous dimensions. The operators with the lowest anomalous dimensions 
are the most relevant in the low-energy sector of the microscopic model, assuming that 
this sector can be mapped onto the free boson theory. 



■^One can easily identify the operators (j>o and with the phase operators, canonically conjugated to 
the number operators N and J. The phase operators are not well-defined in the whole Hilbert space of the 



theory, while their exponents are. Thus, the commutation relations (4-l£) should be understood as applied 
to the operators e""'^" and e™***^ exclusively. 
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We start by considering the following correlation function of the free boson theory 



{{T<j,{x,T)<l,{x' ,t'))) = -^ln[(x - x'f + v\T - r'f]. (4.24) 

Equation (|^) is written assuming that the limit L — > oo is taken. Moreover, to 
remove the ill-defined operator 0o from the fields <j){x,T) and (j){x',T') one should dif- 
ferentiate them at least once with respect to x or r and x' or r'. Thus, an example of 
a well-defined correlation function is 

K (x - x')'^ - v'^(t - t')'^ 



One can see that the correlation function (4.25) exhibits the power-law decay with the 
exponent equal to two as x — x' and r — r' goes to infinity. This exponent defines 
the so-called anomalous dimension (or conformal dimension) of the operator dx(t) ■ it 
is equal to one. Calculating the correlation function {{TdxOdx'O)) one gets that the 
anomalous dimension of the operator dx6 is equal to one as well. The generalization of 
these results is obvious: the anomalous dimension of the operator d'!^d'^(j) is equal to 
n + m. We recall that dxO ~ 9r0j which follows from Eqs. ( |4.7| ) and ([4.^). 

Next, we define two one-parametric families, Vm{x) and Wm(x), of the so-called 
vertex operators, 

Vm{x) =: e^'"'^^^) :, W^(x) =: e*'"''^^) :, (4.26) 
where m is an arbitrary real number. One has 

{{TVrm{x,T)V-m2(.x',T'))) = (4.27) 

{{TWmAx,r)W-m2{x',T'))) = ^-^^ (4.28) 

((rVr„,(x,r)>V 

— r?i2 {x , T 

))) = 0. (4.29) 

Thus, the anomalous dimensions of the operators Vm and Wm are equal to Km? /A 
and m? /AK, respectively. There are no constraints within the free boson theory on 
the possible values of K and m. However, when one uses the free boson theory to 
describe the low-energy sector of a microscopic theory, some constraints can appear. 
In particular, to describe the low-energy sector of the Lieb-Liniger (as well as g-boson) 
model, the possible values of m in ( 4.26| ) should be restricted by the following discrete 
set: 

m = ±l,±2,.... (4.30) 



This condition will become clear from the arguments of sections and p.l| . 
We are interested in the following correlation function of an operator j4(x, r) 

((rA(x, t)A{x' , t'))) as X — x' ^ oo and t — t' ^ oo. (4-31) 

To calculate this function, we assume that the operator A = A{(j)) can be expanded in a 
series and that every term of this series is a (?!)-dependent operator with some conformal 
dimension. The fields with the lowest conformal dimension are either dx(t> and 8-^(1) with 
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the conformal dimension 1 or the vertex operators ( [4.26| ) with the conformal dimensions 
dimVm = Km?/ 4: and dimWm = m?/4:K. Which of these four operators has (or have) 
the lowest dimension, depends on the values of K and m. Thus, we write 

A{(l)) = aid^cj) + a2dr(t> + 03 V^, + a^Wm^ + h.o.t., (4.32) 

where the symbol "h.o.t" stands for the subleading terms, and oi, . . . , 04 are some c- 
number coefficients. We want to stress that the correct usage of the expansion ( [4. 321 ) is 
as follows: among four terms written explicitly on the right hand side, one should select 
the one (ones) with the lowest conformal dimension, and the rest should be included into 
the "h.o.t.". In other words, if, for instance, Vmi has the lowest conformal dimension, 
the first subleading term does not necessarily come from the remaining three terms 
written explicitly on the right hand side of Eq. ( f4.32| ) . 



4.4 Bosonization of the g-boson lattice model: operator 
relations 

In this section we shall exploit an assumption that the low-energy physics of the g^-boson 
lattice model, Eq. ( 2.16| ), is described by the Hamiltonian ( [4.1C1| ). In other words, we 
will exploit an assumption that the g-boson lattice model (and its continuum limit, the 
Lieb-Liniger model) belongs to a universality class usually referred as the Luttinger 
Liquid [12, [l^. This assumption becomes useful in practice when the correspondence 
between the operators of the microscopic theory (the fields Bn and Bn in our case) 
and the operators of the free boson theory is established (shortly, establishing this 
correspondence, one "bosonizes" the microscopic theory). 

Our aim is to bosonize the currents Eq. ( |3.50p . We write them in the form 
given by Eq. ( [4.32| ) and we should find the spectrum of the anomalous dimensions and 
the values of oi, . . . , 04 from the properties of the microscopic theory. To distinguish 
the microscopic currents j*-™^ from their bosonized form, we write the latter as 



Am) 



the current in the microscopic (g-boson) theory 

an operator within the free boson theory corresponding to j 



(i) We associate the operator N appearing in section ^L2 with the particle number 
operator of the microscopic model. The particle number operator commutes with j'-™^ 
Eqs. (|]5|), dl^, and (|3^), while it follows from Eqs. (|1|), (|4l9| ), and ( p^ ) that 
[A^, Wm] 7^ 0. Therefore, the vertex operators Wm are not present in the expansion 
( |432D of the currents 

(ii) We require the bosonized form of j^"^^ to be unchanged under the transformation 
X — !■ x+L. It can be easily seen from Eqs. (4-14) and ( |4.26 ) that for the vertex operators 
Vm this implies the following constraints on the possible values of m : 



m 



0,±2,±4, 



(4.33) 



Upon bosonization, the ground state average (• • • ) should be replaced with the average 
((• • •)). Since (j'-™'^) = 0, one has {(j^^^)) = in the bosonized theory and the operator 



Vm=o = 1 is not present in the expansion ( 4.32| ) of 



•(m) 
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(iii) One can see from (II) and Eq. ( |4.27| ) that the lowest possible anomalous di- 
mension of Vm is K. It can be easily shown, using the results of Ref. |0| , that 



K >1 as 7 > (4.34) 

in the Lieb-Liniger model. For the g-boson lattice model, we shall always assume that 
we are sufficiently close to the continuum limit to satisfy Eq. ( 4.34 ) as well. Taking 
Eq. (4.34) into account, one bosonizes j^"^) as follows 



and 



•(m) 
Jt 



■^^dxcj) + — l9r0 + h.o.t. 
vr vr 



5^0 H 

vr vr 



dx(t> + h.o.t. 



(4.35) 
(4.36) 

(4.37) 



together with the equation of motion ( |4.6| ) imply that the coefficients am and Pm in 
Eq. ( |4.35| ) are the same as in Eq. ( [4.36| ). 
(iv) The coefficients qq and Pq are 



The bosonized form of the continuity equation ( |3.51| ) , 



(m) 



ao = 1, 



and therefore 



;{0) 

Jt 



vr 



/5o = 0, 



-dx(p + h.o.t. 



= -drCp + h.o.t. 

vr 



(4.38) 

(4.39) 
(4.40) 



The proof of Eq. ( [4. 381) will be given in section 5.1 . 

(v) It follows from (iv) that Eq. ( |435D , can be written in terms of the compo- 
nents of the current j ^'^^ : 



Jr 



a. 



+ + h.o.t., m = 0, ±1, ±2, . . . . 



(4.41) 



To use this relation, W6 iiGGci to know the coefficients cxm 

and Pm ■ They will be expressed 
in section 4.5 in terms of the ground state averages of the conserved currents of the 
(7-boson lattice model. 



4.5 Bosonization of the g-boson lattice model: averages 
of the Noether currents 

An important part of the bosonization procedure is the calculation of the non-universal 
c-number parameters of the effective theory (Luttinger parameter K, sound velocity v, 
coefficients Om and Pm) by establishing their correspondence with the properties of the 
microscopic theory under consideration. This is done for the q-boson lattice model in 
the present section. The methodology used is standard for the Bethe-ansatz solvable 



models |14|, so the presentation will be rather brief. 



27 



From the representation (4.41) it is evident that the coefficients and describe 
the response of the system to large scale (smooth) variations of the local density 

■{0} 



and the particle current jx ■ To calculate this response, we take our microscopic model 
and calculate the variation of ^Tr™^ in response to the variation of ^Tr*^^ and J'x^^ in the 
vicinity of the ground state. It follows from Eqs. ( 2.7j ) and (|3.30| ) that 

(JW) = (4.42) 

in the g-boson lattice model (recall that the symbol (• • • ) denotes the ground state 
average). Consider a homogeneous density variation, which satisfies Eq. (4.42). Than 
the coefficient Om in Eq. ( 4.41 ) is defined by the response of {J'r"^^) to this density 
variation: 

(4.43) 



__/ '7(m)\ 



where the density D is given by Eq. (|]7|). To find Pm, we bosonize the operator 
defined by Eq. ^^): 

ie I dx r]{x)j^\x) 



exp 



(4.44) 



and apply to Eq. ( 4.41| ) the same gauge transformation as was discussed in section 3^ 
We thus get 



e=0 



e=0 



+ 3 — 



+ h.o.t.. 



m 



0,±1,±2, 



e=0 



Assuming that r][x) is a linear function of x, 

r/(x) = ax + h, 



(4.45) 



(4.46) 



we use Eqs. ( 3.47 ) and ( 3.49| ) with a playing the role of an infinitesimal variational 
parameter, and finally obtain 



2 



(4.47) 



Using Eq. (2.77) one can rewrite the relations ( |4.43p and ( [4.47 ) in terms of the integrals 
of motion: 

^^'-^ (4.48) 



Otr 



and 



/5„ 



9(/o) 



2 2^ {Im) 



(4.49) 



2 ih) • 

Finally, we calculate the quantity^ Kv using the local gauge transformation gener- 
ated by the operator Q^, Eq. ( [1.44D . Comparing Eqs. ( [4. 7] ) and ( 4.40| ) we find 



-iKvU{x) + h.o.t. 



(4.50) 



In the course of calculating the function (|l.9|) the parameters K and v appear in the combination 



exclusively. 
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This gives the following commutation relation, with the use of Eqs. ([4.3£ 

Kv d Kv 



vr dy 



vr 



-dx5{x - y). 



M) and 



(4.51) 



Now we apply the local gauge transformation generated by the operator Q^, Eq. ( [4.44| ), 
to the x-component of the current j^'^^, Eq. ( 4.40| ). Using Eq. (4.51) we get 



de 



e=0 



de 



i / dyv{y)hfHx),3^^Hy) 



iKv 



IT 



d^rj{x). (4.52) 



Imposing the condition ( 4.46| ) onto r/(x) and comparing Eqs. ( [4.52| ) and ( p. 49 ) we find 
where L, 5, and M are related by Eq. ( p.l9| ). 



5 Contour-independent integral identities 

The bosonization procedure, discussed in section is an effective tool for calculating 
the correlation functions of the microscopic model in the low-energy (long-distance) 
limit. Our task is, however, to calculate the local function ( |1.9| ). In the present section 
we construct some contour-independent integral identities relating the short- and long- 
distance correlation functions of the microscopic model using its integrable structure. 
The long-distance contribution can be found by making use of the bosonization proce- 
dure, and we are thus getting "for free" an infinite set of non-trivial local correlation 
functions. To calculate the function ( |1.9D we need only one equation from this set; 
for other local observables more equations could be needed. It will be also clear from 
the calculations that the method proposed is general and can be used for finding local 
correlation functions within other integrable models. 

In section ^T] we calculate the coefficients ao and /?o, thus proving Eq. ( 4.38| ). 
The method we use in this proof is generalized in section |5.2| to get some nontrivial 



identities relating the conserved currents of the microscopic ((/-boson lattice) model 
with the properties of this model in the bosonized limit. In section |5.3| we use these 
identities to relate some local operator of the g-boson lattice model to the properties 



of this model in the bosonized limit. The result is given by Eq. ( 5.40 ). Its continuum 
limit, studied in section ^, is the exact expression ( 1.12| ) for the function 53(7). 



5.1 Contour-independent integral identities and the bo- 
sonized limit of the density operator 

We consider in the present section a method of calculating the coefficients ao and Pq 



This method provides us with a proof of Eq. (4.38) and, at the same time, it is an 
important component of more general construction studied in section and relating 
the short- and long-distance correlation functions of the microscopic model. 
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Consider the commutator of the vertex operator Wm, Eq. 



with the number 



operator N. This commutator can be calculated using formulas of section 4.2 and gives 



[N,Wmix)] = -mWmix). (5.1) 

The boson annihilation operator iIj{x) of an arbitrary microscopic theory, obeying the 
property 



[iV,V'(x)] = -V(x), 
can, therefore, be written in the low-energy sector as follows 

ip{x) = cWi{x) + h.o.t. 



(5.2) 
(5.3) 



Here c is an unknown constant which depends on the structure of the microscopic 
theory. 

Using Eqs. (p~^), ( |3.13D , and ( 3.50 ), we represent Eq. (5^) as follows: 



dx 



iW(x),V(y)l =-Hy) 



(5.4) 



It would be tempting to substitute the bosonized expressions (4.35) and (5.3) into 
Eq. (|5.4| ) and in such a way obtain an equation for ao and Pq. One should keep in 
mind, however, that the bosonization technique is an approximation working in the 
long-distance limit. Whenever the arguments of two operators, jl-^\x) and ip{y), are 
close to each other, their product ji^^ {x)Tp{y) cannot be bosonized by simply bosonizing 
each of the two operators. To circumvent this difficulty, we use a trick which will play 
a crucial role for our studies of the function (|1.9|). Let us explain this trick in detail. 

To begin with, we recall first the definition of a Green's formula in the classical 
analysis. Consider a domain D in the (x, r) plane (we assume, for simplicity, that V 
is compact and has a piecewise-smooth boundary dV). For any functions P{x,t) and 
Q{x, t) with continuous first derivatives in 2?, the following formula (Green's formula) 
is valid 

: J Pdx + QdT. (5.5) 

V dV 

The boundary dV is oriented counterclockwise: when going along dT> the exterior of 
T) is kept on the right. We now consider the contour integral 



dxdT\ — 

ox OT 



Pdx + Q dr, 



where the contour T is not necessarily closed. If 

dQ _ dP 
dx dr 



(5.6) 



(5.7) 



in some region T>', than, according to Eq. (|5.5| ), any deformation of T within T>' does 
not change the value of the integral ( p.6D . 

As our next step, we perform a sequence of transformations of Eq. (|5.4D . It is clear 
that one can set y = and integrate over the interval (— L/2,L/2) without loss of 
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Figure 1: In (a) we show the integration contour F of Eq. ( |5.8| ). It consists of two horizon- 
tal segments of length L, infinitesimally close to the x-axis. The integration contour F of 
Eq. ( |5.14| ) is shown in (b). It has the shape of a square with side- length L. Both contours 
are oriented counterclockwise. 



generality. Then we write 



fL/2 




fL/2 


/ dx 


= dx 


/-L/2 




J -L/2 



iW(2;,r^ +0)V(0,0) 



V^(0,0)jW(x,r^-0) 



The second argument of the operators j'f^ and is the imaginary time. The symbol 
T denotes imaginary time ordering, defined by Eq. (4.11). The contour T is shown in 
Fig-Sa). 

We now prove that the expression on the right hand side of Eq. ( |5.8D is contour- 
independent. Consider an arbitrary contour enclosing the point (0, 0) in the [x^ vr) 
plane. To make contact with Eqs. ( |5.5D -(^), we write 



dxrj(°)(x,r)^(0,0). (5.8) 



P 



resulting in 



-TjW(x,r)V(0,0), Q = rjW(x,r)^(0,0), 
dQ dP 



dx dr 



a^r#)(^,r)^(o,o). 



(5.9) 



(5.10) 



Recall that j!h^^ is given by Eqs. (|3.15| ) and ( ^.50|) . The T-operator in the expression 



(5.10) commutes with at all points in the {x,vt) plane except at the origin of the 
coordinate system. Therefore 



dQ dP 
dx dr 



ra^jW(x,r)V(0,0) = 0, {x,vt) / (0,0), 
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(5.11) 



where the last equahty is ensured by the continuity equation ( 3.51 ). We have thus 
shown that the expression 



Pdx + QdT 



Tji'>'> (x, r)^(0, 0)dx + Tj"^^^ (x, r)V(0, 0)dT 



(0)/ 



where 



0, 



dx^e^^rj(o)(x,r)V'(0,0), (5.12) 



(5.13) 



is contour- independent, unless the contour crosses the origin of the coordinate system. 
For the contour F plotted in Fig. [l|(a), the equatio n (|5.12| ) reduces to Eq. (p.Sj). We 
have thus proved that the right hand side of Eq. (|5.8D is contour-independent, and 
Eq. d^) can be written as follows: 



(5.14) 



Now we choose the shape of the contour F in Eq. ( ^.14 ) as shown in Fig. ||(b): a 
square with side-length L. For large system size, L, the operators jl^\x, r) and ip{0, 0) 
are well-separated, and their product can be bosonized using Eqs. (g!35|) , (^ ) and 
(|5.3|). Equation (|5.14|) then reduces to a condition 



ao = 1, 



(5.15) 



as already announced in Eq. (4.38). 

The other statement announced in Eq. ( 4.3^ ), /3o = 0, can be proven by considering 
the commutation relation of 'ip{x) with the momentum operator P, Eq. (p..5|): 



We write, using Eqs. ( |3.15| ) and ( p.50D , 



L/2 I- 
-L/2 



if (x),V(o)' 



-- / dxTjf (x,r)V(0,0), 



(5.16) 



(5.17) 



where the contour F is chosen as shown in Fig. |l|(a). Like Eq. (|5.8| ), this expression can 
be written in a contour-independent form, after which the contour can be deformed to 
the shape shown in Fig. |l|(b), and the operators j^^^ and can be bosonized. One gets 
for [P, ^(0)] in the bosonized limit 



-/3ocWi(0)+h.o.t 



(5.18) 



On the other hand, it follows from Eq. ( |5.3| ) that the operator dxi^{x) in the bosonized 
limit is proportional to the operator dx6Wi{x), whose anomalous dimension is higher 
than that of VVi(x). Thus the only way to satisfy Eq. ( ^.16 ) is to require 



/3o = 0. 



(5.19) 



This completes the proof of Eq. ( [4. 38] ) . It is important to stress that the result ( [4.38| ) 
is valid for any interacting system, the only condition that should be fulfilled is the 
existence of well-defined number and momentum operators for the system. 
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5.2 Contour-independent integral identities and the bo- 
sonized limit of the Noether currents 



We have considered in section |5.1| a contour-invariant integral representation of the 
operators [N,'ip{x)] and [P,'ip{x)]. This method estabUshes a connection between the 
short- and long-distance properties of the microscopic theory. The long-distance sector 
of the theory can be bosonized, with the result of getting explicit answers for the 
correlation functions. We will continue to work in the present section with the contour- 
invariant integral representation, studying the Noether currents in the Lieb-Liniger 
model. Recall that we do not distinguish the Lieb-Liniger and g-boson lattice model 
unless the lattice regularization is required explicitly. Therefore, to shorten notation, 
we will use in the most cases the continuum space variable x instead of the discrete 
variable n. The modifications necessary to take into account the discreetness of the 
space variable are obvious. 
We introduce an operator 

W("'\x,t)= r dx' f;^\x',T), (5.20) 

J -L/2 



where is the r-component of the conserved current Eq. ( ^.50 ). The object 
which will play a crucial role in our further calculations is 

Gf'"") = j dx <|ri(")(2;,T)VF(™)(0,0)^ -dr (Tj'i'\x,T)W^'^\<d,Q)^ 

dx^e^^- (ri(")(x,r)Ty(™)(0,0)), (5.21) 



where e^j^u is defined by Eq. ( p. 13 ), and T by Eq. ( 4.11j ). Like in section denote the 



terms in the integrand of Eq. ( ^.21 ) as follows: 

P = - (rjW(x,r)T^("^)(0,0)) , Q = (rjW(x,r)Ty(™)(0,0)) , (5.22) 
therefore, 

^-^ = (9^TiW(x,r)Ty('")(0,0)) . (5.23) 

The T-operator in this expression commutes with dt at all the points of (x, vt) plane 
except the segment 

T = (-L/2 < X < 0,0), (5.24) 

therefore 

^-^ = (ra^jt)(x,r)w(-)(0,0)) = 0, (x,r) ^ T, (5.25) 

where the last equality is ensured by the continuity equation (|3.51| ) . 

Let the contour F in Eq. ( p. 21 ) be a square centered around the origin and with the 
side of length /, as shown in Fig. |a. The segment T of the (x, r) plane is shown there 
by a thick solid line. The cross in Fig. |2|a indicates the intersection point of F and T. 
Any deformation of F not changing the position of the intersection point, leaves the 
value of Gp"'™'^ unchanged. 
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Figure 2: In (a) we show the integration contour T used in section |5]^. It has the shape of a 
square with side-length L. The integration contour V used in section p.3| is exhibited in (b). 



We want to replace the exact currents j^"*^ entering Eq. (5.21) by their approximate 
expressions (4.35) and ( 4. 361 ) obtained within the free boson theory. There is, however, 
a problem: the operator jl^\x,T) and the operator ji"^\x',0) entering W^"^\0,0) 
are not separated by the asymptotically large space-time interval in a vicinity of the 
intersection point of F and T, and one cannot apply Eqs. ( 4.35| ) and ( 4.36 ) to the 
operator j^"'^(x,r)jr"*^(x',0). To circumvent this problem we use the following trick: 
Introduce an operator 

/•L/2 

W^"'\x,t) = - dx' j^;^\x',T). (5.26) 



Combining Eqs. ( |5.20| ) and ( [5.26| ), one gets 



L/2 
-L/2 



dx'j|™)(x',- 



{In 



(5.27) 



tion 2.3. Then, we split the contour F into two parts: 



where Im are the integrals of motion of the g-boson lattice model, discussed in sec- 

(5.28) 



F+ X > 
F_ X < 



Using Eqs. ( 5.26D -( ^.28 ) we rewrite the expression (5.21) in the following form 



G 



(n,m) 



dx^e^^" (Tjt^(x,r)T^("^)(0,0) 



+ 1 dx^e^-' (Tji^\x,T)W^^\0,0)) + j dx^e^^''(Tji^\x,T){Im-{Irn))). 

(5.29) 
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The third term on the right hand side of Eq. ( ^.29| ) vanishes because the ground state 
|gs) is an eigenfunction of Im- In the first (second) term the distance between the 



local field {x, t) and the local field jr {x' , 0) entering the operator (0, 0) (the 
operator ^^^"^^(0,0)) is larger than /. 

A nontrivial output from Eq. ( ^.29| ) comes in the limit 



oo, 



oo, 



/ 

L 



0. 



(5.30) 



Since the operators jl^\x,T) and jr"^\x',0) entering Eq. (|5.29| ) are separated by a 



distance larger than /, one can use the bosonized expressions (4.35) and ( [4.36| ) in the 
limit ( ^.30 ). Taking into account Eq. ( 4.24| ), one gets after some algebra 



G 



{n,m) 



Kv 



TT 



(5.31) 



5.3 Contour-independent integral identities: the main re- 
sult 

In present section we relate the ground-state average of some nontrivial local operator 
of the g-boson lattice model to the properties of this model in the bosonized limit. Es- 
tablishing this relation, Eq. ( |5.4[1| ), we get all necessary information for the calculation 
of the function 93(7), Eq. (R 



We use the techniques developed in sections 5J and 5^. Choose the contour V in 
(x, vt) plane as shown in Fig. ^(b). Namely, V consists of the two horizontal segments, 
(— / < X < e, +0) and (— / < x < e, —0), and the square with the side /, centered around 
the origin. The latter part is just the contour T used in sections 5.1 and 5.2. Using the 
techniques developed there we write 



dx 



i(")(x,0),W^(™)(0,0) 



0. 



(5.32) 



The parameter e > plays the role of a regularization parameter, as it is necessary 
to treat unambiguously the behavior of the integral in the vicinity of x = 0. We shall 
work with Eq. ( |5.32| ) in the limit ( |5.30| ), and we can therefore use Eq. ( 5.31 ) to get 



dx 



dx' 



L/2 



f^Hx),i^\x' 



Kv 



vr 



(5.33) 



We remind the reader that we write the lattice regularization explicitly only when 
it is necessary. Here is the proper place to do it. Equation ( ^.33 ) then takes the form 



k=-l5~^ p=-L5-^/2 



Kv 



vr 



Let us take the limit ( 5.30 ), keeping 6^^ finite. Therefore 

L5^^ — > 00, 16~^ 00 
in Eq. ( 5.34| ). We then choose the parameter e such that 

e6^^ > m. 



(5.34) 

(5.35) 
(5.36) 
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We set n = —1 and m = 2 in Eq. (5.34). The local operators ji (k) and j^^ (p) are 
defined by Eqs. ( 2.45| ), ( |2.46| ), ( |2.51| ), and (p.50| ). Using these definitions, we calculate 
the expression on the left hand side of Eq. ( |5.34| ) under the conditions ( 5.35| ) and ( |5.36| ): 



^ ^ ^ fc=-M-ip=-L5-i/2 



A Rt nt , 



Comparing this result with the right hand side of Eq. ( 5.34 ) one gets 



(5.37) 



+ x^b]b]b,B,^^) = -1(1-^) ^(a_i/32 + a2/3_i). (5.38) 



r 

2 



Kv 
vr 



This expression relates the ground-state average of the local operator of the microscopic 
model to the properties of this model in the bosonized limit. 
Let us combine Eqs. (|5]3|), (|1|), and ( ^) 



{b]b]b,b,^,) = + ^ (i - y) v^"-^^^ + 

and substitute into Eq. ( ^.39|) the results ( [4.48|) , (|4.49| ) and (|4.53|) . This gives 

5 (A) 



2; xHW (/i)5(/o) 



(5.39) 



(5.40) 



6 Three-body local correlation function 

Working with the g-boson lattice model, we obtained in previous sections all the for- 
mulas necessary to calculate 53(7). The remaining task is to take the continuum limit 



in these formulas and to collect all them together. The limit is taken in section |6.1| and 
the formulas are collected together in section |6^. 



6.1 Properties of close to the continuum limit 

In this section we continue our studies of (7^) in the g-boson lattice model, started in 
section 2^. The continuum limit ( |2.19| ) of the g-boson lattice model gives the Lieb- 
Liniger model, Eq. ( |1.1D . The limit in which we are interested is the continuum limit 
together with the extra condition L — > 00, coming from Eq. ( |5.30| ). Therefore, we 
consider 

5^0, M^oo, K^O, L^oo, (6.1) 



where 



L = 5M, 



c/2 = Kb 



-1 



(6.2) 



Let us renormalize the quasi-momenta pj used in section 2.4 



(6.3) 
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The kernel ( |2.86| ) can be represented in the hmit (6.1) as follows 
K{p) = 6-^K{k), K{k) 



When written in terms of the variables kj, the normalization condition (2.84) becomes 

r.A(5) 



n 



A(5) 



dk p{k, 5), 



and the Lieb-Liniger equation ( ^.85 ) becomes 



Pik,5)-- 



dkK{k - k)p{k,5) = —. 

-A{5) 



27r 



(6.5) 



(6.6) 



Here the ground-state density D is defined by Eq. (2.79). Note that the functions A 
and p in Eqs. ( |6.5D and ( |6.6D are different from those used in Eqs. ( 2.84 ) and ( |2.86| ). 
We, however, use the same symbols for these two couples, since we shall work with 
Eqs. (^) and (|6.6D exclusively. 

The ground state expectation values of /i and I2, Eqs. ( 2.87 ) and ( p. 88 ), can be 
represented in the limit (|6.1D as follows 



L 



xHD--h2 + -h4 + ... 



and 



where 



L 



2! 



4! 



hm{6) 



dkp{k,5)k'" 



(6.7) 

(6.8) 
(6.9) 



Equations (6.7) and ( |6.8D are not series expansions in powers of 6 since the functions 
hm themselves depend on 5 through the quasi-momentum distribution p and through 
A. The series expansion of hm in powers of 5 is 



Let us consider Eqs. ( |6.5D and ( |6.6| ) at 5 = : 

r-A(O) 



and 



p(fc,o)- — 



D 



1 rm 



A(0) 



dkp{k,0) 



dk- 



2c 



-Pik,0) = — . 
2^7-A(0) c'^ + {k-k)^^^ ' 27r 



Following Ref. ||^] we change the variables 

A; = A(0)z, c = A(0)a, p(A(0)z, 0) = cr(z). 



(6.10) 



(6.11) 



(6.12) 



(6.13) 
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When written in these variables, Eqs. (6.11) and (6.12) become Eqs. ( |1.15| ) and (1.14), 
respectively. The function h!^ introduced by ( |6.10 ) can be written as 



htHD,c)=D^+^em{l). 



(6.14) 



where the function €^(7), given by Eq. (|1.13| ), depends on the dimensionless parameter 
7, Eq. (p.. 7]), only. One verifies readily that satisfies the following differential 
equation 

+ ^^ht^ - ("^ + 1)^^^^ = 0. (6.15) 



Next, we consider the solution of Eqs. ( |6.5D and ( |6.6| ) to second order in 5. Substi- 
tuting Eq. (6^) into (|6.6|) and keeping the terms up to the order of 5"^ we find 



1 rm 



dk- 



2c 



1 



pik,6) = -(l + —D 



2vr J-A(5) + {k- ky 27r 
By rescaling the quasi-momentum distribution function 

c5^ 



c5^ 



p{k,6) = ^l + —D]p{k,6) 
we find that p{k, 6) satisfies the integral equation 



(6.16) 



(6.17) 



p{k,S)-- 



1 fMS) 



dk- 



2c 



^^J^m c^ + {k-kY^ ' 27r 



(6.18) 



with the condition ( |6.5| ) renormalized as follows: 



-m 



dk p{k, 6). 



(6.19) 



It is clear from comparison of Eqs. ( |6.19| ) and ( |6.18| ) with Eqs. ( |6.11| ) and ( 6.12| ) that 
the function /i2, Eq. ( |6.9| ), can be represented to order 5^ as follows 



h2{D,c) = 1 + 



6 



A(5) 



dkk'p{k,S) = 1 + 



6 



h'^' D 



6 



Expanding this equation and dropping terms of higher order than 5^ , we get 



6 



(6.20) 



(6.21) 



Finally, combining Eqs. ( |6.15| ) and ( 6.21| ), we get the following differential equation for 
dh2 dh2 , cS'^D f dh2\ , , 



where we have dropped all the terms of the order higher than 5"^ 
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6.2 Three-body local correlation function: the result 



We are now in the position to derive the main result of this paper: equation ( [1.12 ) for 
the three-body local correlation function (73(7), Eq. (|1.9| ), of the Lieb-Liniger model 
in the limit (1.11). To do this, we take the limit (6.1) in the ground-state averages 
of local operators of the g-boson model, which were obtained in earlier sections. The 
corresponding calculations are straightforward but rather lengthy, and we shall only 
sketch their main steps. 

The desired ground state expectation value (1.9) in the Lieb-Liniger model can be 
obtained by taking the limit ( |6.lD of the following g-boson local field: 



(1 - x')B\B\B,B,+r +(1-^1 <li-.{B]B,+^ 



d 



dq 



c6^ 
3 



+ 



(6.23) 



where the dots denote terms of higher order in 6. To prove Eq. ( |6.23D we use Eq. ( |2.96| ) 
and then take the limit ( |6.1| ) in the way discussed in section |2.2| , Eqs. ( p.22 )-( |2.25| ). 
After rather lengthy algebra we got the right hand side of E q. (|6.23| ) . 

As our next step we rewrite the left hand side of Eq. ( |6.23 ) using the identities 



(5.40) and (O 



1 - x')b]b]b,B,^, + ( 1 - ? ) l4-iB]B,^i 



2 

1+(1 



dq' J 
2 



X^ 



1 d 
1 — ] Q — 



d{h) 
d{h) 



2 J^dq\ L 
{h)d{Io 



(6.24) 



and take the limit (|6.ip in the resulting expression. Substituting the expansions (pT 



and ( p. 81) into the right hand side of Eq. ( 6.24 ) one arrives at 



l-x'')BjBjBjB,+,+ {l 



2 



dq 



(B B 



53 f dh2 , dh2 



3h-, 



12 r 5^+'^^ 



3h 



1 

Vie 



15/l4 + 



h2 dh2 



(6.25) 



where we have dropped terms of higher order than S'^. Using Eq. ( |6.22| ) to transform 
the first term on the right hand side of Eq. ( 6.25| ) one can see that the expansion 
( 6.25 ) starts from the terms of order 5^. Comparing these leading order terms with the 
term written explicitly on the right hand side of Eq. ( 6.23] ) and replacing hm with 
Eq. dill), one gets after some algebra the final result, Eq. ( |1.12| ). 
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